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The form factors for the semileptonic decays of heavy-light pseudoscalar mesons of 
the type D Kev are studied in quenched lattice QCD at /3 = 6.0 using Wilson fermions. 
We explore new numerical techniques for improving the signal and study 0(a) corrections 
using three different lattice transcriptions of the vector current. We present a detailed 
discussion of the relation of these lattice currents to the continuum vector current and 
show that the disagreement between the previous results is to a large extent due to the 
value of Zv used in the calculations. We also present results for the decay constants of 
light-light, heavy-light and heavy-heavy mesons. 
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1. Introduction 

The physics of mesons containing one heavy valence quark {c. h) and one Hght valence 
quark (w. d. s) is of considerable interest at present. Experimental investigations of these 
systems may provide the most accurate determination of the parameters of the standard 
model responsible for flavor mixing and CP violation (the CKM matrix).* This can only 
be achieved, however, if we are able to ol)tain a quantitative understanding of the influence 
of the strong interaction on their structure and decays. 

Numerical simulations of lattice QCD provide a solution to this problem, and there 
has been considerable activity in the field in recent years. Two groups have presented 
results for the form factors for the semi-leptonic decays of heavy-light pseudoscalar and 
vector mesons, namely Bernard et al. [2] [3] [4], and the Rome-Southampton group [5] [6] 
[7]. These results have large statistical errors, and, in certain instances, are in conflict. To 
resolve these discrepancies and check for systematic errors, further studies using different 
methodologies are necessary, and this is what we undertake in this paper. 

In this paper we present data on the pseudoscalar and vector meson decay constants 
for light-light, heavy-light and heavy-heavy systems. We then move on to the main topic 
of this paper semileptonic form factors for pseudoscalar meson decays such as Z) — >^ Kev. 
The structure of this paper is the following. We review the phenomenology of semi-leptonic 
form-factors in Section 2. In Section 3 we discuss the renormalization of lattice operators 
and give details of the lattice setup in Section 4. Results for the decay constants are 
presented in Section 5 along with those for the vector current renormalization Zy- The 
analysis of the form-factors is given in Section 6 along with a discussion of the quality of 
the signal and a comparison with previous data is given in Section 7. Finally, we state our 
conclusions in Section 8. 

2. Phenomenological Background 

We consider the case. D^Xlv. where X has flavor content us {K or K*). In the one 
W exchange approximation the amplitude is 



A recent overview of the experimental status of charmed meson physics is given in ref. [1]. 
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where Gp is the Fermi constant, Vcs is the c s CKM matrix element. This process 
is particularly simple because the hadronic and leptonic currents factorize. The leptonic 
part of the decay can be calculated accurately using perturbation theory, while to take 
into account non-perturbative contributions to the hadronic part 

= (X|^7^(l-75)c|i^) (2.2) 

one resorts to lattice QCD. In this paper we present results for the simpler of the two 
cases, i.e. — )• K~e^v. 

The matrix element iJ^ can be parameterized in terms of two form factors: 

{K-{pK)\s^f,{l - 75)c \D\pn)) = pMQ^) + (lf^f-{Q^)^ (2-3) 

where p = (pj^ + p^) and q = {pp) — px) is the momentum carried away l)y the Icptons, 
and = —q^ (which is always positive). We use the Euclidean notation p = (p,iE) so 
that p^ = p"^ — E"^ . An alternative parameterization is 

(Jf-(pjr)|37„(l-75)c|D°(pD)) 



where 



/o(g') = f+m + o o f-{Q^). (2.5) 



In the center of mass coordinate system for the lepton pair, i.e. g = or equivalently 
Pk = Pd , one has 

{K-{pK)\s^c\D\pn)) = 2pnf+iQ^), 

{K-{pKp^,c\D\pn)) = ^^^^^MQ'). ^^'^^ 

Thus, the form factor /+(Q^) is associated with the exchange of a vector particle, while 
/o(Q^) is associated with a scalar exchange. It is common to assume nearest pole domi- 
nance and make the hypothesis 



f+iO^) — ~. — {^i/^ 2 ' /o(Q^) — — {^lo 2 ' (2-''') 



where nijp is the mass of the lightest resonance with the right quantum numbers to mediate 
the transition; Z)+(1969) or Z)*"'"(2110) in the pseudoscalar or vector channels respectively. 



The goal of the lattice calculations is to determine the normalizations /+(0) and /o(0) and 
map out the dependence. 

In the limit of vanishing lepton masses, the vector channel dominates and one can 
write the the differential decay rate as 



Gl\v... 

1927r3m35 " ' ' ' " - ■ " ■ ^2.8) 
MQ^) = (ml + ml- Q^ f - Am%m\. 

To integrate this, the functional form of must be known. Assuming vector meson 
dominance, numerical integration gives 

T{D^ K-e+u) = 1.53|V;.n/+(0)|^ x IQ-^sec-^ (2.9) 

Using Eqn. (2.9) one can extract Vc^ once T(D^ — > K~c'^p) has been measured and 
/+ calculated using lattice QCD. In this case D*^ — > K~c~^p. however, \Vcs\ = 0.975 is 
known very accurately, so one extracts |/+(0)| ~ 0.75. The quantity /o(0) has not been 
determined. 

The present study, whose goal is to investigate different numerical techniques in order 
to improve the signal to noise ratio, shows that first principles calculation of form-factors 
can be carried out reliably with today's massively parallel computers. 



3. Renormalization of lattice operators 

Lattice transcriptions of continuum operators like the vector and axial currents suffer 
from 0(a) corrections. To get a handle on these 0(a) effects in the study of form factors, 
we use three transcriptions for the vector current 

V';'^^Hx) = q,(x)lMx), 

V^^'^'ix) = i {qi(x)-f^U^(x)q2(x + a^) + qi(x + afi)-ff,Uf,(x)^ q2(x)) , (3,1) 
y^^°"^-(aj) = i {q^(x)(-ff, - l)Uf,(x)q2(x + a^i) + q^(x + ai^)(-ff, + l)Uf,(x)^ q2(x)) . 

We remind the reader that V^'^^^'(x) is conserved only for degenerate quarks. In the calcula- 
tion of the pseudoscalar decay constant we use only the local axial current ^^(35)757^^2(35). 
In these currents the quarks qi and q2 may both be light, or one heavy and one light, or 
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both heavy. The first step in relating these lattice currents to their continuum counterparts 
is to calculate the normalization of lattice fields. 

Lepage and Mackenzie have given a prescription for relating lattice quantities to their 
continuum counterparts [8]. They show how to handle the large tadpole contributions 
to operators and as an example show that the renormalization constant Zy for the local 
vector current agrees very well with non-perturbative measurements even for ma k, 1 
[9]. Since the reliablity of calculations of matrix elements depends crucially on accurate 
determination of these parameters, we present an explicit 1-loop analysis and reproduce 
the results of Lepage and Mackenzie. 

The overall philosophy of the 0{a,^) improvement scheme is that in perturbation 
theory tadpole contril)utions are large when quantities are expressed in terms of the l)are 
parameters. The three l)asic ways in which tadpoles contril)ute to Feynman diagrams are 
shown in Fig. 1. These can and should be absorbed into the renormalized quantities, 
so that the resiilting perturbation theory is well behaved. The goal is to absorb the 
correction to the fermion line in the definition of the renormalized mass, to the gluon line 
in the definition of the renormalized coupling, and whenever possible to cancel the tadpole 
contributions in operators by combining them with the renormalized quantities. We now 
show how this works to order for the different currents we use in this paper. 

3.1. Field renormalization 

We want to calculate Z^, the relative normalization of lattice fields, defined as 



On dimensional grounds, Z can only be a function of r, ag and m.a for Wilson fermions, 
where ma is the qiiark mass in lattice units. The inverse of the Wilson propagator is (we 
set r = 1 to simplify the discussion) 



where the factor of 1/a has been extracted to make it dimensionless. The 1-loop con- 
tribution / has the structure / = Sq -|- Si^a -|- S2ma and it contains the tadpole con- 
tribution (i$a — 4)K [10] [11]. The Sj and K are numerical constants. In this equation 
ma is defined to be the bare value 1/2k — 4. We shall keep track of terms of 0{ag) and 
0(ma), but neglect 0(a1) and O(asma) terms. Thus, the results are strictly speaking 
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valid only for agina < a^, though in practice the range of validity has to be determined 
non-pert urbatively. 

We now set A; = and express quantities in terms of = ik^. For brevity of notation 
we write f(—ikQ, ma) as f{kQ, ma). Then the forward moving part of the inverse propagator 
is 

{2kV) ^ a = — expkoa + 1 + ma + asf{koa,ma). (3-4) 

The renormalized mass is defined by the location of the pole in the propagator, which is 
given by the relation 

ex.p k^a = 1 + ma + agf (kQa,ma). (3-5) 
The zero mass limit is determined from the condition {mc = 1/2kc — 4) 

mca + asf{0,mca) = (3-6) 

which is used to define the value of Kc 

1 



aj{0,mca) 



c 



2k^ 

8kc = 1 + ^asf{0,mca) 



(3.7) 



Close to the chiral limit ma is of order ag as can be seen from Eq. (3.6). In order to 
calculate the wave-function renormalization one needs to evaluate the residue at the pole. 
This is 



, df{koa,0) , n\ df{koa,0) 

expkoa - — = 1 + ma + a;^/(fco«,0) - — 

okna okna 



[1 -|- ma + asf{kQa, 0)][1 — a. 



Of 

df{koa,0) . 
dkna 



(3.8) 



where we have set ma = in terms proportional to as these are 0{a^^ma). Note that 
all the terms proportional to ag have contributions from tadpole diagrams. In lattice 
calculations it has been traditional to factor the residue into two terms, the first term is 
a lattice effect and defined to be the relative renormalization of the lattice to continuum 
field ip and the second is the square of the wave-function renormalization that has to be 
included for each external fermion leg in the operator. The seemingly bad behavior of 
perturbation theory arises if these two terms are truncated at different orders of ag. We 
consider it useful to give explicit details of how the tadpole contributions are distributed 



and absorbed into the renormalized quantities so that the resuhing perturbation expansion 
is weh-behaved. 

The first step in the reorganization is to renormahze the bare quark mass and express 
it in terms of k^. Using Eq. (3.7) we get 

1 



1 + ma + agf{kQa, 0) = 3 + agf{kQa, 0) 

2k 



AKr 



3k 

AKr 



+ as fikoa,0) 



8k 



(3.9) 



Now the term proportional to is of order a^ma, and can therefore be dropped. This 
rearrangement gives the final result for the residue of the propagator V (note that we have 
to multiply Eq. (3.8) by 2k), i.e. 



8Kr 



3k 

AKr 



d 



a. 



' dkna 



f{ha,0) 



(3.10) 



which is to be evaluated at the pole. The tadpole contributions are distributed as follows. 
As shown in Eq.(3.7) the factor 8kc absorbs 1/4 of the tadpole term that is not proportional 
to ^ . This is equal and opposite in sign to the tadpole term proportional to $ that is 
included in the standard wave-function renormalization (^l — asdf{kQa, 0) /dkooj . Therefore 
these two pieces always cancel. We reiterate this observation in slightly different words to 
emphasize the point. The tadpole contribution to 8kc in the perturbative expansion 



= 1 + 1.36Aav + ... 



(3.11) 



is equal and opposite to that from wave-function renormalization of the external legs. This 
is true for any operator O. Thus the combination 8kcZo will bo better behaved at 0{as). 



The remaining tadpole part is included in ^/l — 3Ki/AKc. On using the measured value for 
Kc this factor becomes independent of perturbation theory and is therefore well behaved. 
This completes our discussion of tadpoles on the fermion line. A way to include tadpoles 
on gluon lines will be discussed later when we describe the Lepage-Mackenzie scheme. 

Using this scheme we can now write down the normalizations for the vector current. 
The 1-loop calculations of the local vector and axial currents have been done by Martinelli 
and Zhang [10] and we have checked their calculations. The local vector current contains 
a tadpole contribution from the wave-function renormalization, so we combine it with the 
perturbative expansion for 8kc to get 



qi{x)jf,q2{x) 



cont 



cont 



1 - 



3ki 

AKr 



1 - 



3k2 

AKr 



(1 - Q.82av)qi{x)-i^^q2{x) 



(3.12) 



Here ay is the renormalized coupling that is defined later. To get the normalization for 
extended currents one has also to include the 0(ma) terms in the vertex. For the conserved 
current the vertex at tree level is modified to e~^°"7p. The extra factor e~^o" has been 
obtained by making the same simplifications as in Eq. (3.4). This factor exactly cancels 
the [1 + ma + asf{kQa, 0)] part of Eq. (3.8). The remainder (1 — Si) is equal and opposite 
in sign to the 1-loop corrections to the vertex. Thus the renormalization factor for the 
"conserved" current is unchanged from the tree level definitionf 



cont 



(3.13) 



The above result also applies to the case when the current is fiavor changing but its matrix 
elements are taken in the forward direction. On the other hand the vertex for the extended 
1-link current is 7^ cos k^^a, which reduces to 7^ at the order to which we are working. The 
two tadpole contributions, one from the wave-function renormalization 1 — Si and the 
other from the expansion of the link in the operator are again equal and opposite. Since 
these two cancel one docs not need to combine the perturbation result for Z^^*' with that 
for 8Kr. Thus we use 



cont 



= - ^Jl - ^(1 - im8ay)V^^'- 



AK.r V 4k^ 



(3.14) 

L 



It should be emphasized that our analysis does not specify whether it is better to use the 
perturbative or non-perturbative value for Sk^ in cases where the tadpoles in the pertur- 
bative part cancel. The difference between using the perturbative and non-perturbative 
residt is fti 6% at (3 = 6.0. i.e. SKf^^*' = 1.185 versus the non-perturbative value of 1.256. 
This difference is a measure of the residual uncertainties in the renormalization of the 
currents. 

Lastly, the vertex for the local axial current gets a tadpole contribution only from 
(1 — Si). This can be cancelled by combining it with the perturbative expansion for Sk^. 
Then 

(3.15) 



cont » ^'^c V ^f^c 



Lepage and Mackenzie have given a mean-field prescription for determining the re- 
lation between lattice and continuum quantities [8]. The only difference in the final re- 
sults for local operators in the mean-field approach is to replace the prefactor 8kc by 



f We thank G. Martinelli for discussions on this point. 
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where 



^piaq is the trace of the plaquette. However, within the mean- 



field framework SkcC/q = 1 and deviations from unity are a measure of uncertainty in this 
approach. To give the reader a feel for how consistent these approaches are we give the 
numerical values at /3 = 6.0. They are I/C/q = 1.140 and '^kJJq = 1.1. The perturbative 
resuh for I/Uq is 1 + lM7av = 1.165. 

To summarize the Lepage-Mackenzie prescription for removing the potentially large 
tadpole contributions in local fermionic operators at 0(a,) the cook-book recipe is 
(a) For each quark of flavor i in the operator use the normalization 



between the lattice field and its continuum counterpart, 
(b) Calculate the corrections to the operator O. l)oth on the lattice and in the contin- 
uum in say the MS scheme. In the result Zq replace a^are by ay where [8] 



This replacement removes the large tadpole contril)utions to gluon lines that would 
otherwise show up at O(a^) and make the series look badly behaved. 

(c) Determine the scale q at which ay is defined (one prescription for doing this is given 
in Ref. [8]) and scale o;y(7r/a) to avig) rising 2-loop running. With these three things 
in hand one has the final relation Zq between the lattice and continuum operator 
to 0{as). In the rest of the paper we shall choose ay = 5'jj/47r as the renormalized 
coupling at scale q = lAn/a, where g'j^ = l-Tgl^re at /? = 6.0. 

(d) To extend the analysis to non-local operators one has to include the 0{ma) tree level 
corrections to the operator. These have to be combined with the renormalization of 
the field ijj at the stage of Eq.(3.8). This can change the relation between the lattice 
and continuum field ijj from that given in Eq. (3.16) as we have shown by the examples 
of the two 1-link vector currents. 

The lesson from our analysis is that one cannot apply the mean-field prescription to 
calculate the renormalization factors for extended currents by simply counting the number 
of ijj fields and links in the operator. Instead one has to work self-consistently to a given 
order in ag and include 0{ma) and 0{pa) terms present in the Taylor expansion of the 
operator. 



cont 






(3.17) 
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We shall present results using the improved normalizations given above and make 
comparisons with the naive scheme (without the improvement incorporating tadpoles on 
fermion lines), i.e. 4'lont — V2k* ■i/'^. In both cases, however, we will use the boosted g"^ in 
ay to sum up the tadpole contribution to the ghion lino. To allow the reader to evaluate 
the difference between the two normalizations we give the normalization factor for both 
schemes and their ratio in Table 1. As the numbers show, the difference due to the choice 
of the normalization of the currents grows rapidly with the quark mass. 

Our perturbative analysis is only valid close to the chiral limit and takes into account 
0(ma) effects only for on-shell quarks and for processes with no momentum transfer, i.e. 
it is not 0(a) improved. These conditions are not realized in the semi-leptonic decays 
of heavy quarks. Since the normalization factors have a large uncertainty especially at 
large cjuark mass, it is necessary to check the perturbative renormalization factors by 
non-perturbative methods where possible. El Khadra et al. [9] have measured Zy^"'' by 
calculating the matrix element {J/'i/j\Vl"'^'"''' \J/ip) for heavy quarks (m^a ^ 1) and find 
that Eq. (3.12) works incredibly well. This encouraging result needs to be verified for a 
variety of different matrix elements and we discuss this issue further in Section 5. 

3.2. Extension to Sheikholeslami-Wohlert Action 

The above analysis can be easily extended to improved actions. The case of the 0{a) 
improved action proposed by Sheikholeslami and Wohlert [12] is particularly simple. This 
is because the tadpole contributions to the propagator with the SW action are the same 
as those for the Wilson action. Secondly, the perturbative expansion for is unchanged. 
Thus the expression (3.10) is still valid except that the extra fermion-gluon vertex changes 
the function / and consequently the precise form of 1 — Si. Similarly, the finite parts of 
the 1-loop result for the vertex are also different. These corrections have been calculated 
in Ref. [13] and all one needs to do is extract the appropriate perturbative corrections to a 
given operator from it. The rest of the renormalization factors are the same as for Wilson 
fermions and are given by the steps (a,c,d) descril)ed aliove. 

The numerical data show that the improved normalizations work better for the SW ac- 
tion as a consequence of 0(a) improvement. For example, at /3 = 6.0 the non-perturbative 
value for the chiral limit has been measured to be Kc = 0.14556(6) [14]. Thus 8kc = 1.1645 
agrees very well with the perturbative result 1.185 and SkcUq = 1.02. 

9 



4. Details of the lattice setup 

4.1. Lattice parameters 

Our statistical sample consists of 35 lattices of size 16^ X 40 at /3 = 6.0. This sets of 
lattices have been used previously for spectrum and weak matrix element analysis using 
both Wilson and staggered fermions [15] [16] [17], and the details of the lattice generation 
are given in [15]. 

In this study our main goal is to investigate different numerical techniques in order 
to improve the signal to noise ratio. For this purpose we use only one value of the heavy 
quark mass, k = 0.135, and only two values of the light quark mass, k = 0.154 and 

0. 155. Our results for meson masses expressed in lattice units are given in Table 2. Using 

= 1.9 GeV, these corresponds to a heavy meson of mass 1.54 and 1.59 GeV (about 
the mass of the physical charm quark) and to light-light pseudoscalar masses of roughly 
690 MeV and 560 MeV. Our heavy-light pseudoscalar mesons therefore correspond most 
closely to the physical D meson, with a somewhat massive light constituent, while the light- 
light mesons are analogous to the physical K. We will henceforth adopt this nomenclature. 

4.2. Quark propagators 

The calculation of quark propagators is done on lattices doubled in the time direction. 

1. e. 16"^ X 40 — > 16"^ x 80. We use periodic boundary conditions in all four directions. 
These propagators on doubled lattices are identical to a linear combination of propagators 
calculated with periodic (P) and antiperiodic (A) boundary conditions on the original 
16^ X 40 lattice. For the source on time slice 1. the forward moving solution (time slices 
2-40) is F = (P + A) /2 while the backward moving solution (time slices 80-42) is B = 
{P-A)/2. 

The details of our implementation of the "Wuppertal" smeared source method [18] 
are described in [16] and we have reused the light-quark propagators generated in earlier 
calculations. For a given light quark mass we have four types of quark propagator which 
we denote G^^(x; 0), G^g{x; 0), G^gg{x; 0) and G'^g(x; 0) for "light local-smeared", "heavy 
local-smeared", "light smeared-smeared" and "heavy smeared-smeared" . The "smeared- 
smeared" quark propagators are obtained by applying the Wuppertal smearing procedure 
at each sink time slice of the "local-smeared" propagators. 

For the calculation of 3-point functions we require a fifth type of quark propagator, 
namely a light-heavy propagator with the insertion of a zero 3-momentum smeared pseu- 
doscalar source at some fixed time Iq. This propagator, which we denote G^j, is calculated 
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by doing a heavy quark inversion using a zero-momentum solution of the hght propagator 
on a specified time-shce as the source. In terms of quark propagators without insertions 
we have the definition 



G'i^{x-o) = Y,Glsi^-^y^toh5G'ssiy,to-,o)- (4.1) 



Note that wo perform "Wupportal" smearing twice at to in the construction of G^'. This 
is needed to cancel the smearing factors when constructing ratios of 3-point to 2-point 
correlators. We choose to = 32 in order to be far from the boundary at t = 40. Thus 
wrap-around effects in time direction are exponentially damped by at least 18 time slices 
because we use propagators calculated on doubled lattices. 

Our method for extracting form-factors is similar to that used by the Rome- 
Southampton group [5] [6] [7] except that we have used smeared quark propagators which 
improve the signal and calculate the form-factors for three different currents. In Sections 
6 and 7 we comment on the improvements over previous results and make a detailed 
comparison. 

4.3. 2-point and 3-point correlators 

The large time behavior of a 2-point correlator at a given 3-momentum is 

X 

^ ^2WM^^^ Ol |0), t ^ oc (4.2) 



exp -tJmfo +P 



where Oi and O2 arc interpolating operators for the hadronic state, \h) of mass Mh 
{Ehipf" = p"^ -\- M"^), and Zi is the amplitude for creating and annihilating the state. The 
subscript i in Zi stands for both the flavor of the quarks and the type of smearing for the 
source and sink. We have collected data using both the local and smeared pseudoscalar 
density and axial vector current 

P{x) = qi{x)'y5q2{x), 

(4.3) 

Af^{x) = qi{x)-ff^-f5q2{x), 
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to interpolate for pseudoscalar mesons. The three point function giving the matrix element 
of current, J^, between hadronic states is defined to be: 

(4.4) 

- ^MnMP) ^2(0))(/i2(0) J^\h,{p)){hi{p)\0\\Q), 

for < t < to- 

In order to extract i?^ we construct the following two ratios of 3-point to 2-point 
correlation functions. For the case Oi = 02=P and = V^''^"^ in Eqs. (4.2) (4.4), and 
using the LS heavy- light and light-light 2-point correlators, the ratio TZls{p) is defined to 
be 

Q"^(p,t) 

Tlrqiv) = 

C UP,t)C^LUO^t o - 1) (4 5^ 

- V ^LS (P) Zl's{^){K{p) I I7m c\d{^)), 0<t<to. 
Similarly, using the SS heavy-light and light-light 2-point correlators gives 



T^ss{p) 



^ssiP^ ^)^ss(^^ ^0 ~ ^) (4.6) 



4s(p)^ss(0) {K{p)\s^^c 



D{0}), 0<t<t 



0- 



The meson decay amplitudes in Eqs. (4.5) and (4.6) are determined in separate fits. The 
results of these fits are combined at the time of jackknife analysis to extract i?^. This 
procedure is described in more detail in the next sub-section. Note that the two ways 
of calculating differ only in the 2-point correlators used to cancel the exponential 
dependence on t and subsequently on which meson decay amplitudes have to be removed. 
Thiis the difference lietween the results is a measure of the imcertainty in the fits. We 
show results for the two cases separately and for the average of the two as it reduces the 
fluctuations. The normalization of the currents is incorporated at the very end of the 
analysis as it is a common overall factor. 

4.4. Fitting procedure 

To extract amplitudes and masses from the large time behavior of 2- point correlators 
we assume that at large t only the lowest mass state dominates the correlation functions. 
To ensure this we first examine the effective mass plot for a plateau and then make a single 
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mass fit selecting the range of the fit based on the following criteria: (a) t^^^ always lies 
in the plateau, (b) tmax is selected to be as large as possible consistent with a signal. 

In most cases we find that the central value obtained from the fits is the same with 
and without using the full covariance matrix. In some cases we cannot use the full range 
of the plateau because the covariance matrix is close to being singular. In such cases the 
problem is not that we cannot invert the covariance matrix but that the resiilt is very 
sensitive to the range of the fit (the central value can change by one or more standard 
deviations on the addition of a single point to the fit range) and/or some of the jackknife 
samples give grossly different results. Our tests using suljsets of the data and by examining 
the covariance matrix show that this instability is a result of inadequate statistics. Some 
of the causes for this bad behavior are discussed in Ref. [19]. To overcome this problem we 
reduce the range of the fit until we get a stable result and compare it with that obtained 
without the full covariance matrix to check for consistency. 

All error estimates are obtained using a single-elimination jackknife procedure. Pre- 
vious analysis [20] leads us to believe that the lattices used are sufficiently decorrelated for 
this method to be adequate. 

In order to extract the matrix element iJ^ we make two separate fits for each jackknife 
sample to (1) the ratios of correlators defined in Eqs. (4.5) and (4.6) assuming that the 
residt is constant over times slices on which the lowest state satmates the 2-point corre- 
lation functions; (2) the 2-point correlators needed to extract the extra decay amplitudes 
that have to be removed from TZls and 7^55. The results of the fits are combined for each 
jackknife sample and the final value and error on it is calculated by the jackknife method. 
This process is explained in more detail in Sections 5 and 6. 

5. Meson decay constants 
5.1. Pseudoscalar mesons 

The definition of the pseudo-scalar decay constant on the lattice is [21] 



where Z'^Za is the axial current renormalization, and /-^ = 132 MeV in our convention. 
The details of the four different methods we use to extract the decay amplitude are given 



ZIZa{Q\AY' 



Tr{p)) 
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in Ref. [16], so we only present the results here. All four methods give consistent estimates 
of /tt and we quote their mean in Table 3 for three different choices of Z'^Za- 

Since we combine results from different correlators in order to extract /^r we select the 
fit range based on the following criteria: (1) goodness of the fit, (2) presence of a plateau 
with a similar mass estimate from each of the individual correlators. To monitor for the 
possibility that the covariance matrix is not well determined we make fits with and without 
the full covariance matrix. In cases where the results with the covarience matrix deviate 
from the data we decrease the range of the fit to sec if the results stabilize, otherwise 
we quote the central value from the uncorrelated fits and give the larger of the two error 
estimates. 

Results for /-^ at p = (0, 0, 27r/16) are also given in Table 3 and are consistent with 
p = data. The qiiality of the signal is good enough to extract reliable numbers even 
though the plateaus in the effective mass are much shorter than those for the -p = (0, 0, 0) 
case. At higher momenta the errors get much larger and the determination of /-^ is not 
reliable. 

Some of the data in Table 3 have been published previously in Ref. [16]. The old 
numbers may differ from those given in Table 3 under label C because we reanalyzed 
all the data and the subjective choice of the range of the fits may be different. Thus 
the difference between the two results should be treated as a measure of the associated 
systematic error. 

5.2. Vector mesons 

The vector decay constant is defined in terms of the matrix element of the local vector 
current Vi between the rho and vacuum states 



hi our calculation V- is used to both create and annihilate the vector meson. The two 
methods used to extract fy^ are described in Ref. [16] and we quote the mean value of 
the two values in Table 3 using three different choices for Z3,Zv. 



(o|y/-- \p) 




ZlZvfv 
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5.3. Analysis 

The data in Table 3 show that estimates for the decay constant have a very large de- 
pendence on the normalization of the currents. Results for with the Lepage-Mackenzie 
scheme deviate by a large amount from experimental values as shown in Fig. 2 (case A in 
Table 3). Choosing a different value for the at which ay is to be evahiated changes the 
numbers by < 5%. which is insignificant compared to the difference at the heavier masses. 
The results using the non-perturbative value of Zy arfi in much better agreement (case C 
in Table 3 and Fig. 2). The non-perturl)ative values of Zy. given in Tal)le 4. are calculated 
using the ratios of 2-point correlators for rho mesons, for example Zy'^' is defined as 

m = o)vi"-{p = o,t))' ^"-'^ 

with Z^ = V2K* for V^^"'^^'. The data for Zy^*' is obtained using a similar ratio with V/"'^' 
replaced by Vf'^^'. Note that in the calculation of using the above non-perturbative 
value of Zy, the factors of Z^ in cancel, so one cannot test a particular normalization 
scheme by this procedure. 

In the case of the zero momentum correlator the vector current does not insert any 
momenta at the sink time-slice, so there should be no 0{a) corrections to the conserved 
current. This is why Eq. (5.3) is considered a good way to define Zy. The data in Table 4 
has been obtained iising Z^ = V2k* for Vl"'^'- and TA^'^*- to allow comparison with previous 
pul)lished work. We find that Zy'^' and Zy'^" show very little dependence on the quark 
mass. Furthermore, the results are essentially imchanged when Zy is calculated from 
correlators at non-zero momentum. To convert to a different normalization, say Lepage- 
Mackenzie, recjuires multiplying the numbers in Table 4 by the ratio of the respective Z^ 
factors, in which case Zy develop a dependence on the mass. 

The same problem arises in the calculation of form-factors. The perturbative and 
non-pertmbative normalizations give significantly different resiilts and it is the improved 
perturbative normalization that seems to give consistent results for semi-leptonic form 
factors as we will disciiss below. This dependence of the normalization on the nature of 
the initial and final state in the matrix element has been exposed and investigated by the 
Rome-Southhampton group [22]. In a recent study of semi-leptonic form factors they [23] 
have argued that this mismatch is a symptom of 0(a) effects and have presented data 
showing that the differences decrease on using the Sheikholeslami-Wohlert 0(a) improved 
action. Clearly this issue needs further study in view of the fact that the differences are 
so large. 
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6. Semileptonic form factors for D K e+i/ 

To orient the reader we first summarize the ways in which our method is technically 
different from those used in previous calculations. 

1. We use smeared quark propagators to improve the signal while the previous calcula- 
tions used point sources. 

2. In our approach the interpolating operator for the final state meson sits at the source 
point (t = 0) of the original smeared light propagators. The source for the initial D 
meson is at t = 31 and obtained by making a smeared pseudoscalar insertion in the 
light-to-heavy propagator G^"^^ . By using quark propagators calculated on doubled 
lattices we essentially eliminate wrap-around effects. With the two interpolating op- 
erators at a fixed separation, we can carry out the final contraction for the matrix 
element by placing the weak current at all intermediate times < t < 31. At each 
time slice t we average the matrix element over all spatial points in a state of definite 
momentum transfer. One advantage of this method is that non-local lattice currents 
can be used to calculate the matrix elements. This approach is same as that used by 
the Rome-Southhampton group [5] but different from that used by Bernard et al. [3]. 

3. We calculate the matrix elements for three different transcriptions of the lattice current 
on the same set of lattices. Bernard ci al. have only used v'""''^ while the Rome- 
Southhampton collaboration used only y^"" '*-. 

Our calculation should be regarded more an exploration of methods rather than pro- 
viding hard numbers that can be compared to experimental data because of the number of 
approximations made. We find that the data is of much better quality for the case where 
the final meson is annihilated by the pseudoscalar density ■075 as compared to the axial 
current V'7475V'- For this reason we only present results for the pseudoscalar case. 

In figs. 3-5 we show data for the ratios TZ^s '^ikI T^ss for the three lattice transcriptions 
of the vector current. We have multiplied the IZss data by 10^ and TZ^s by —1 to put them 
on the same plot. The plots show that the signal is very similar with all three currents. 
When making fits to the data we find that TZss has a better plateau and fits can be made 
using the full covariance matrix. The same is not true of TZls ^nd these results have been 
obtained without including the correlations. 

To extract a composite quantity like the form-factor we have combined results from 
separate fits for each jackknife sample as our data sample is too small to make a simul- 
taneous fit to all the necessary correlators. For example, to calculate we first make 
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separate fits to the ratio TZss for V4 and Vi and to the 2-point correlators Cgg and Cgg 
needed for the decay ampHtudes occuring in Eqs. (4.6). A second estimate is similarly 
obtained by using Eq. (4.5). The difference between the two is a measure of the statistical 
uncertainty in the data. A third estimate is constructed by averaging the TZls ^i^d TZss 
results for each jackknife sample. 

The results for all three estimates are given in Tables 5-7 using the normalizations 
defined in Eqs. (3.12), (3.14) and (3.13). There is a ^ 20% difference between the results 
obtained from TZss ^nd TZls for the case of momentum transfer p = (7r/8,0,0). This 
is a manifestation of the fact that the signal in non-zero momentum correlators falls very 
rapidly with p and this presents the largest source of error in the calculation of form-factors. 
We regard the third estimate (average of LS and SS) as our best estimate, however, the 
reader should be aware that there is no deep reason for this. 

We find that all three lattice transcriptions of the vector current give consistent results. 
The results with are about 10% higher. This difference is of the same magnitude as 

that between the perturbative and non-perturbative value for 8kc and would be removed if 
we use the perturbative value for it in Eq. (3.14). It is this consistency which suggests that 
it is the improved normalizations for the currents that work well for the 3-point matrix 
elements. 

The results for matrix elements in the forward direction are, within errors, 1.0 for all 
three currents. For non-zero momentum transfer our result for /_|_ is consistent with the 
phenomenological estimate /_|_ « 0.75, however, we have used too few values of quark mass 
and momentimi transfer to determine whether pole-dominance is an accurate description 
of the variation of /_|_ with respect to Q^. We also find that /_|_ « /o while the errors in 
/_ are too large to make useful deductions. 

If we express results using the normalization that seems to work for /p, i.e. the im- 
proved with the non-perturbative value for Zy, then the results are no longer as 
consistent. For example the results for /_|_ in Table 7 for ki = 0.154 get changed to 0.37(7), 
0.42(11) and 0.60(9) for V/^"'^-. V/^'^^- and V^""-"- respectively. The difference in results with 
the two schemes is siibstantial and indicates that the 0(a) effects with Wilson fermions 
are large and that the improved perturbative normalizations may not be very reliable for 
all matrix elements at /3 = 6.0. 
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7. Comparison with previous results 

There are two previous caleulations against which we can make direct comparisons 
as these use almost the same lattice parameters. The group of Bernard ct al. [3] have 
measured the form-factors on 24^ X 40 lattices at the same values of (3 and k. They used only 
the local vector current, and adopted yet a different normalization. Converting their result 
to our "improved" normalization gives fo{p = 0) = 0.85(10) at k = 0.154 to be compared 
with our value of 0.91(9). Similarly the Rome-Southampton group [6] [7] have measured 
the form-factors on 20 x lO'^ x 40 lattices at the same value of (3 and similar k. They use 
the "conserved" vector current with the same normalization for the vector current that we 
use. Interpolating their results to k = 0.154 and noting that their momentum (q = 27r/10) 
(and hence Q^) is slightly different we get /_|_ = 0.67(5) and /o = 0.65(4), to be compared 
with our results 0.60(9) and 0.62(10) respectively at (q = 27r/8). 

Thus we find that all three calculations give consistent results. The central remaining 
issue is to understand how to reduce the large uncertainty in the normalization of the 
currents. 

8. Conclusions 

We have given explicit details of how to self-consistently derive renormalization factors 
for fermionic operators to 0{as). Our analysis reproduces the Lepage-Mackenzie mean- 
field results for local operators. We also show how to extend the analysis to non-local 
operators and to the Sheikholeslami-Wohlert improved action. 

The results for decay constants and form-factors show large differences depending on 
the choice of normalization used for the currents. No one scheme seems to give reliable 
results for all the different matrix elements. These differences grow rapidly with the quark 
mass. Further work, possibly at weaker coupling and/or using 0(a) improved actions, is 
required to resolve this issue. 

We find an acceptable signal in the 3-point correlators needed to calculate semi- 
leptonic form-factors iising quark propagators with smeared sources. We show that differ- 
ences in previous results is a consequence of different renormalization constants used for 
the vector current. Using the improved normalizations, we find that all three transcriptions 
of the vector current give consistent results. The biggest limitation in the calculation of 
form-factors is the rapidly decreasing signal in the non-zero correlators. Based on present 
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data we believe that major improvements in results for the form-factors are possible by a 
high statistics study on larger lattices, say 32^ x 64. 
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Ki = 0.155 


K,i = 0.154 


Ki = 0.155 


= 0.154 

ny2 — l^.XOO 


Ki = 0.135 


V Kb 1 1 LjJ . j 


231 


235 


270 


272 


316 


V^^^inaive) 

V / 


0.218 


0.217 


0.204 


0.203 


0.190 


Ratio 


1.058 


1.084 


1.326 


1.342 


1.661 




280 


286 


328 


331 


384 


V / 


0.266 


0.265 


0.249 


0.248 


0.232 


Ratio 


1.052 


1.079 


1.319 


1.335 


1.653 




310 


308 

v7 . (J\70 


289 


288 


270 




0.310 


0.308 


0.289 


0.288 


0.270 


Ratio 


1.000 


1.000 


1.000 


1.000 


1.000 


A^^^iimp.) 


0.249 


0.253 


0.291 


0.294 


0.340 


A^°'^(naive) 


0.240 


0.238 


0.224 


0.223 


0.209 


Ratio 


1.039 


1.065 


1.302 


1.318 


1.631 



Table 1: The normalization constants for the three lattice transcriptions of the vector 
current and for the axial current used in this study for the 5 different comV)inations of 
fjuark masses. We give the improved and "naive" normalizations along with the ratio of 
the two for ease of comparison. 
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Ki = 0.155 
i^r, — n 1 fifi 

A.2 — U.idO 


Ki = 0.154 
(•„ — n 1 fizt 

K2 — U.i04: 


Ki = 0.155 

fi 2 — U.iod 


Ki = 0.154 

fv2 — U.iOd 


Ki = 0.135 
— n 1 

K2 — U.iOd 


TT^^ (0,0,0) 
TT^S (0,0,0) 

Trf (0,0,0) 

TT^S (0,0,0) 


0.299(9) 
0.300(8) 
0.297(16) 
0.303(13) 


0.362(6) 
0.365(7) 
0.369(10) 
0.370(7) 


0.813(5) 
0.812(4) 
0.815(8) 
0.812(6) 


0.833(5) 
0.833(4) 
0.836(7) 
0.832(6) 


1.216(4) 
1.217(3) 
1.217(4) 
1.216(4) 


7rSS(l,0,0) 
7r^S(l,0,0) 
Trf (1.0.0) 
7rP(1.0.0) 


0.5(1) 
0.5(1) 
0.5(1) 
0.5(1) 


0.49(3) 
0.51(3) 
0.52(3) 
0.51(3) 


0.89(1) 
0.89(1) 
0.90(1) 
0.90(1) 


0.91(1) 
0.91(1) 
0.92(1) 
0.92(1) 


1.265(7) 
1.267(4) 
1.268(7) 
1.267(5) 


(0,0,0) 
(0,0,0) 
pSS(l,0.0) 
P^^ (1,0,0) 


0.407(16) 
0.411(10) 
0.6(1) 
0.6(1) 


0.459(14) 
0.460(8) 
0.60(3) 
0.62(2) 


0.842(8) 
0.840(6) 
0.93(2) 
0.93(2) 


0.862(7) 
0.862(6) 
0.94(2) 
0.94(1) 


1.229(4) 
1.230(3) 
1.278(8) 
1.280(5) 


af (0,0,0) 
af (0,0,0) 
af (1,0,0) 
af (1,0,0) 


0.7(1) 
0.7(1) 
0.7(1) 
0.8(1) 


0.7(1) 
0.7(1) 
0.8(1) 
0.8(1) 


1.08(5) 
1.11(4) 
1.2(1) 
1.2(1) 


1.07(4) 
1.10(3) 
1.1(1) 
1.2(1) 


1.42(3) 
1.44(2) 
1.47(3) 
1.49(3) 



Table 2: Meson energies extracted from 2-point correlators. LS and SS denote local- 
smeared and smeared-smeared corrletors respectively. Two types of correlator were used 
for the to extract the pseudoscalar mass: (a) (PP) denoted tt: and (b) {A4A4) denoted 
7T2- The energy is given for the particle at rest (the mass) and for the case of one unit of 
momentum. Quoted errors are statistical. 
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Ki = 0.155 
K,2 = 0.155 


Ki = 0.154 
K,2 = 0.154 


Ki = 0.155 
K2 = 0.135 


Ki = 0.154 
K2 = 0.135 


K-i = 0.135 
K2 = 0.135 


(A) A(0,0,0) 
(A) /^(1,0,0) 

V / J /IV" " / 

(A) f -Ho. 0. 0) 


0.076(5) 
0.084(6) 
0.34(2) 


0.089(5) 
0.090(5) 
0.327(17) 


0.125(4) 
0.126(5) 
0.195(6) 


0.131(4) 
0.133(5) 
0.200(6) 


0.206(4) 
0.209(6) 
0.201(4) 


(B) ^(0,0,0) 
(B) f^(l,0,0) 
(B) /-^ (0,0,0) 


0.073(7) 
0.080(7) 
0.33(2) 


0.084(6) 
0.085(5) 
0.31(1) 


0.098(4) 
0.098(5) 

0.148(5) 


0.101(3) 
0.101(4) 
0.150(5) 


0.126(4) 
0.129(5) 
0.121(4) 


(C) /.(0,0.0) 
(C) /.(1, 0,0) 

(c)/-i(o,o,o) 


0.082(7) 
0.090(7) 
0.27(2) 


0.093(6) 
0.094(5) 
0.25(1) 


0.107(4) 
0.108(5) 
0.119(4) 


0.111(3) 
0.112(4) 
0.121(4) 


0.141(5) 
0.143(5) 
0.098(4) 



Table 3: The meson decay constants fp and fy extracted from LS and 5' 5' 2-point corre- 
lators. Quoted errors are statistical. The three cases differ in the normalization used for 
the currents. Case A uses the improved normalizations defined in this paper; case B uses 
i^cont — V2K* and perturbative values of Za and Zy with boosted g"^; case C uses 
'^cont ~ V2k* ip'^j^, Za = 0.86 and the non-perturbative estimate Zy = 0.57. 
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Ki = 0.155 
= 0.155 


= 0.154 
K2 = 0.154 


K-i = 0.155 
K2 = 0.135 


Ki = 0.154 
K2 = 0.135 


Ki = 0.135 
K.2 = 0.135 


Z|7-(p = 0) 

V V-r II y 

z|r-(p=i) 


0.572(5) 
0.55(2) 
0.58(2) 


0.571(3) 
0.57(1) 
0.56(1) 


0.565(2) 
0.544(8) 
0.562(6) 


0.565(3) 
0.545(7) 
0.563(5) 


0.551(2) 
0.536(5) 
0.548(3) 


Z;^>^t-(p*= 0) 

^r-(p=ll) 


0.685(5) 
0.68(2) 
0.69(2) 


0.686(3) 
0.68(1) 
0.70(2) 


0.695(2) 
0.692(5) 
0.690(7) 


0.696(3) 
0.693(4) 
0.693(6) 


0.706(2) 
0.705(3) 
0.704(3) 



Table 4: The renormalization constants Zy' and Zy' extracted from ratios of LS 2- 
point correlators for the p meson. These numbers have been obtained using Z^ = \^2k 
normalization for all three vector currents. We give the results for zero momentum, and 
for one unit of momentum l)oth parallel and perpendicular to the index of the current. 
Quoted errors are statistical. 
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K = 0.154 


Current 


/+(g2 = -0.043) 


/_(g2 = -0.043) 


/n(g2 = -0.043) 


/o(g2 = 0.222) 


'(r Local 


0.56(13) 


-0.33(23) 


0.58(14) 


0.91(12) 


-{/■Ext. 


0.60(12) 


-0.34(24) 


0.62(14) 


0.98(11) 


'(/'Cons . 
V 


0.55(11) 


-0.20(16) 


0.56(11) 


0.88(13) 


K = 0.155 


Current 


/+(g2 = -0.021) 


/_(g2 = -0.021) 


/o(g2 = -0.021) 


/o(g2 = 0.268) 


'(/'Local 


0.59(19) 


-0.54(30) 


0.61(20) 


0.91(9) 


'(/'Ext. 


0.60(23) 


-0.43(31) 


0.61(24) 


1.00(10) 


V 


0.55(20) 


-0.22(25) 


0.56(21) 


0.87(9) 



Table 5: Form factors for the semi-leptonic decay D — )■ K~e'^p using the ratio defined 
in Eq. (4.6). The results for /_|_, /_ and /o in columns 2-4 are for momentum transfer 
q = (7r/8,0,0) while column 5 gives /q at § = 0. Quoted errors are statistical. 
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K = 0.154 


Current 


/+(Q2 = -0.043) 


/_(Q2 = -0.043) 


/n(Q2 = -0.043) 


/o(g2 = 0.222) 


'(r Local 


0.67(12) 


-0.56(31) 


0.73(15) 


0.91(10) 


-{/■Ext. 


0.76(18) 


-0.47(30) 


0.81(18) 


1.02(15) 


'(/'Cons . 
V 


0.65(10) 


-0.25(20) 


0.68(11) 


0.89(9) 


K = 0.155 


Current 


/+(g2 = -0.021) 


/_(g2 = -0.021) 


/o(Q2 = -0.021) 


/o(g2 = 0.268) 


'(/'Local 


0.71(23) 


-0.76(42) 


0.78(24) 


1.01(13) 


'(/'Ext. 


0.73(25) 


-0.61(41) 


0.78(25) 


1.10(14) 


V 


0.67(22) 


-0.33(33) 


0.70(21) 


0.98(12) 



Table 6: Form factors for the semi-leptonic decay D ^ K e'^v using the ratio defined in 
Eq. (4.5). The rest is as in Table 5. 
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K = 0.154 


Current 


/+(g2 = -0.043) 


/_(g2 = -0.043) 


/n(g2 = -0.043) 


/o(g2 = 0.222) 


'(r Local 


0.61(11) 


-0.45(25) 


0.66(13) 


0.91(9) 


-{/■Ext. 


0.68(13) 


-0.41(24) 


0.72(14) 


1.00(12) 


'(/'Cons . 
V 


0.60(9) 


-0.23(17) 


0.62(10) 


0.88(9) 


K = 0.155 


Current 


/+(g2 = -0.021) 


/_(g2 = -0.021) 


/o(g2 = -0.021) 


/o(g2 = 0.268) 


'(/'Local 


0.65(20) 


-0.65(36) 


0.69(21) 


0.96(10) 


'(/'Ext. 


0.66(24) 


-0.52(36) 


0.70(24) 


1.05(11) 


V 


0.61(21) 


-0.28(29) 


0.63(20) 


0.93(10) 



Table 7: Form factors for the semi-leptonic decay D ^ K e'^v using the average of 
and SS data. The rest is as in Table 5. 
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Fig. 1: The tadpole contributions to (a) fermion line, (b) gluon line, and (c) operator. 



28 



0.4 



0.3 - 



I > 



0.2 - 



0.1 




0.2 



0.4 0.6 



0.8 1 



Fig. 2: The vector decay constant as a function of m^/m^. The experimental numbers 
are labeled by the symbol □. The data given as case A (C) in Table 3 are marked by the 
symbol x (o). 
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Fig. 3: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the local vector current for (a) V4{p = (0,0,0), (b) Vi{p= (7r/8,0,0)), and 
(c)y4(p = (vr/8,0,0). 
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Fig. 3: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-slice using the local vector current for (a) V4{p = (0,0,0), (b) Vi{p= (7r/8,0,0)), and 
(c)y4(p = (vr/8,0,0). 
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Fig. 3: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the local vector current for (a) V4{p = (0,0,0), (b) Vi{p= (7r/8,0,0)), and 
(c)y4(p = (vr/8,0,0). 
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Fig. 4: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the extended vector current for (a) V4{p = (0,0,0), (b) Vi{p = (7r/8,0,0)), 
and (c) V4(p = (7r/8,0,0). 
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Fig. 4: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the extended vector current for (a) V4{p = (0,0,0), (b) Vi{p = (7r/8,0,0)), 
and (c) V4(p = (7r/8,0,0). 
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Fig. 4: The data and fit to TL^s and TLss versus the seperation from the kaon's annihilation 
time-shce using the extended vector current for (a) V4(p = (0,0,0), (b) Vi{p = (7r/8,0,0)), 
and (c) V4(p = (7r/8,0,0). 
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Fig. 5: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the conserved vector current for (a) V4{p = (0,0,0), (b) Vi{p = (7r/8,0,0)), 
and (c) V4(p = (7r/8,0,0). 
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Fig. 5: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the conserved vector current for (a) V4{p = (0,0,0), (b) Vi{p = (7r/8,0,0)), 
and (c) V4(p = (7r/8,0,0). 
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Fig. 5: The data and fit to TZls ^^id TZss versus the seperation from the kaon's annihilation 
time-shce using the conserved vector current for (a) V4{p = (0,0,0), (b) Vi{p = (7r/8,0,0)), 
and (c) V4(p = (7r/8,0,0). 
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